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Abstract 

The Pachmarhi Array of Cerenkov Telescopes consists of a distributed array of 25 
telescopes that are used to sample the atmospheric Cerenkov Photon showers. Each 
telescope consists of 7 parabolic mirrors each viewed by a single photo-multiplier 
tube. Reconstruction of photon showers are carried out using fast timing information 
on the arrival of pulses at each PMT. The shower front is fitted to a plane and the 
direction of arrival of primary particle initiating the shower is obtained. The error in 
the determination of the arrival direction of the primary has been estimated using 
the split array method. It is found to be ~ 2.4' for primaries of energy > 3 TeV. 
The dependence of the angular resolution on the separation between the telescopes 
and the number of detectors are also obtained from the data. 

Key words: VHE 7 - rays, Atmospheric Cerenkov Technique, angular resolution, 
non-imaging telescope array, PACT. 



1 Introduction 



The 7— ray sky has been successfully observed by satellite based missions in 
the energy range between 50 keV and 30 GeV [1,2]. At higher energies the 
steeply falling 7— ray spectrum makes it impossible to observe by satellite 
based detectors, primarily because of the limited detector size and exposure 
time available for a given source. Alternatively, the ground based atmospheric 
Cerenkov technique has been successfully exploited at these energies. Over 
the years, this technique has proved to be the most sensitive in studying the 
celestial 7— rays in the very high energy (VHE) range, ~ 100 GeV - 50 TeV. 

Celestial VHE 7— rays initiate an electromagnetic cascade in the atmosphere as 
they enter the earth's atmosphere. The electrons and positrons in the cascade, 
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being relativistic, emit Cerenkov light as they propagate down the atmosphere 
resulting in a faint flash of bluish light lasting a few nanoseconds. This direc- 
tional Cerenkov flash is detected on the ground by conventional optical light 
detectors during moon-less clear nights. Taking advantage of the fact that the 
Cerenkov flash is highly collimated, within a cone of half angle ~ 1°, along 
the direction of the incident particle most experiments simply limit the optical 
field of view of the Cerenkov telescopes to a small region of the sky. 

However the main drawback of this technique is the presence of abundant 
cosmic ray background which severely limit the sensitivity of this technique. 
Numerous efforts to develop methods to distinguish the Cerenkov light pool 
produced by cosmic 7— rays from that by the cosmic rays led to two impor- 
tant techniques based on complementary ways of viewing the cascade, viz. the 
angular sampling or the imaging technique and spatial sampling or the wave- 
front sampling technique. Both these techniques are currently being employed 
by different groups [3-6]. 

The possibility of timing the Cerenkov light front for the estimation of arrival 
direction of the primary has been realized as early as 1968 by Tornabene 
and Cusimano [7]. It was later demonstrated by the Durham group [8] that 
while viewing Crab pulsar the on-axis events were richer in 7-rays than off- 
axis events. It was also used by Gupta et al.[9] to improve the signal to noise 
ratio of atmospheric Cerenkov telescopes consisting of large mirrors with poor 
optical quality. The basic principle of all these experiments is the fast timing 
technique using spatially separated Cerenkov telescopes. 

The signal to noise ratio in such an experiment is given by [10] 

S I AT ^ 
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where F c and F 1 are the fluxes of omni-directional cosmic rays and 7— rays 
from a point source respectively. A is the collection area of the array, T is the 
time of observations, uj is the solid angle of acceptance, f r is the fraction of 
showers due to cosmic rays that are identified and rejected as background and 
/ 7 is the fraction of showers due to 7— rays that are identified as signal and 
hence retained. In order to achieve high -j| apart from increasing / 7 and f r one 
could either increase the collection area and the observation time or decrease 
the solid angle of acceptance. For a given exposure time and available hardware 
resources one can possibly increase S/Nby only reducing uj as celestial 7— rays 
from point sources are directional while cosmic rays are isotropic. Due to the 
finite opening angle of the Cerenkov cone and the spread in the arrival angle 
of Cerenkov photons the aperture of the telescopes have to be restricted to few 
degrees, which sets a lower limit to uj. However, it is possible to improve the 
4r for point sources without losing Cerenkov light if the direction of arrival of 
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primary particles is estimated accurately [7,9,11,12]. The shower axis retains 
the original direction of the primary. 

To estimate the arrival direction accurately, the error in the reconstructed 
direction, i.e. the angular resolution has to be very small. The two dominant 
factors which contribute to the angular resolution are the average distance, D, 
between the telescopes and St, the uncertainty in the measurement of arrival 
time of photons at the telescopes. For example, the angular resolution in the 
zenith angle, 9 for n detectors is given by [9] : 

= (2) 

Dcos9\n K J 

Therefore, a large number of telescopes with instrumentation to measure the 
relative arrival time of photons and separated by large distances are needed 
to reconstruct the shower front and estimate the direction of arrival of the 
shower fairly accurately. On the other hand, in case of the imaging technique 
the angular resolution is limited by PMT (pixel) size which is typically of the 
order of a quarter degree. The improvement in the signal to noise ratio by 
restricting the angle of acceptance to S9 would be: 




This is a very significant advantage of non-imaging arrays with high angular 
resolution. Moreover, this method of reducing background is independent of 
the primary species. Hence it is extremely useful at very low primary energies 
where the cosmic ray electrons form a significant source of background. Elec- 
trons cannot be discriminated easily, unlike hadrons, since they too undergo 
electromagnetic interactions in the atmosphere akin to 7— rays. However one 
has to bear in mind that the angular resolution could be poorer at lower 
primary energies due to paucity of Cerenkov photons at the ground level. 

In the rest of the paper we discuss the method of obtaining the arrival direction 
of the primary radiation using our array of Cerenkov telescopes and estimate 
its error. 



2 Pachmarhi Array of Cerenkov Telescopes 

The experiment at Pachmarhi (longitude: 78° 26' E, latitude: 22° 28' N and 
altitude: 1075 m), is based on the wavefront sampling technique and employs 
an array of Cerenkov telescopes, called the Pachmarhi Array of Cerenkov 
Telescopes(PACT), to sample the Cerenkov light pool. PACT is now fully 
functional. 
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The experimental set-up of PACT has been explained in detail elsewhere [13]. 
Briefly, it consists of a 5 x 5 array of atmospheric Cerenkov telescopes deployed 
in the form of a rectangular matrix with a separation of 25 m in the N-S 
direction and 20 m in the E-W direction. Each telescope consists of 7 parabolic 
mirrors of 0.9 m diameter mounted paraxially and having a focal length of 90 
cm. Each mirror is viewed by a fast photo-multiplier tube (PMT, EMI 9807B) 
at the focus behind a circular mask of ~ 3° diameter. However the field of 
view is limited by the diameter of the PMT photo-cathode to ~ 2.9° FWHM. 

The array has been divided into 4 sectors with six telescopes in each. The 
pulses from 7 PMTs in a telescope are added linearly to form a telescope sum 
pulse called royal sum. Each royal sum from the 6 telescopes in a sector are 
suitably discriminated (typical royal sum rates ~ 30-50 kHz.) and a trigger is 
generated by a coincidence of any 4 of these 6 royal sums. The typical event 
rate is ~ 2-5 Hz. For every event, information on the relative arrival times and 
density of Cerenkov photons are recorded for the 6 peripheral mirrors/PMT 
in each telescope in each sector. The relative arrival times of royal sum pulses 
are recorded both in the respective sector and in the central data processing 
station. 

Thus, PACT measures the arrival time of shower front at various locations 
within the Cerenkov light pool at two distance scales, short range (intra- 
telescope) and long range (inter-telescope) . The arrival direction of the shower 
is estimated from the measured arrival times of Cerenkov photon front at each 
of the spatially separated telescopes while the six adjacent measurements in a 
given telescope could be used to study the fluctuations in the measured shower 
parameters. The dispersion of photon arrival times at each telescope contain 
the signature of the primary species [14] and hence could be used to distin- 
guish between 7— rays from the background [15]. The density measurements 
on the other hand enable us to estimate the energy of the primary species as 
well as to reject cosmic ray background [16]. 



3 Alignment of Mirrors and Telescope 
3.1 Pointing Accuracy of Telescopes 

The telescopes are equatorially mounted and each telescope is independently 
steerable in right ascension and declination within ±45°. The movement of 
the telescopes is remotely controlled by Automatic Computerized Telescope 
Orientation System (ACTOS) [17]. The hardware for ACTOS consists of a 
semi-intelligent closed loop stepper motor drive system which senses the an- 
gular position using a gravity based angle transducer called clinometer with 
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an accuracy of V . The two clinometers, one each in N-S and E-W direction, 
are accurately calibrated using stars at various hour-angles and declinations. 

In order to estimate the source pointing error of the system, the telescopes 
were oriented to different bright stars at random positions in the sky and the 
offsets in orthogonal directions (ascension and declination) with respect to the 
star were noted. Using this data it is concluded that the system can orient 
the telescopes to the putative source with a mean offset of 0°.003 ± 0°.l. The 
source pointing is monitored with an accuracy of 0°.05 and corrected in real 
time during tracking. The uncertainty in pointing represents the subjective 
error in manual reading of the star position at the cross wire of the guiding 
telescope field of view (± 0°.05). 



3.2 Measurement of Alignment Accuracy of Mirrors 

As the seven mirrors of a telescope are mounted on a single mount it is neces- 
sary to ensure that all their optic axes are parallel to each other so that they 
view the same part of the sky. This alignment of mirrors is done manually by 
positioning the star image at the centre of the focal plane. 

To check the accuracy of alignment of the mirrors and telescopes a bright 
star drift scan is carried out. The telescope is aligned to an isolated bright 
star (typically of visual magnitude 2 to 3). Then the telescope is moved to 
the west by 4° and at a suitable time the telescope tracking is switched off. 
The counting rates from each of the PMTs are monitored every second and 
recorded. The count rates stay reasonably constant until the star walks into 
the field of view when they increase, go through a maximum and return to 
the background value as shown in figure 1. The background count rates before 
and after the star transit are fitted to a linear function. The background 
is subtracted from the count rates during the star transit by interpolation. 
The background corrected count rates are then fitted to a quadratic function. 
Figure 1 shows one such fit to a typical count rate profile of a mirror. Table 1 
shows the summary of results of a typical bright star scan. It shows in column 
3 the FWHM of the drift scan profile of the count rates due to the star. The 
offset of the centroid of this profile with respect to the centre of the field of 
view {i.e. expected transit time) is shown in the last column. This offset is 
the error in the alignment in right ascension. The relative FWHM's of the 
profiles could indicate the misalignment in declination, if any. On the other 
hand the absolute values of FWHM of the count rate profiles are a function of 
the PMT gains, image quality (point source image size < 1°), star brightness 
etc. Using this method, it is ensured that the optical axes of all the 7 mirrors 
in a telescope are parallel to each other within an error of 0°.2 or less. If the 
error exceeds this value the particular mirror is re-aligned and re-checked. 



5 



2500 



2000 



o 1500 



o 
O 



1000 - 



500 



17.0 



X 



/ 



/ 



/ 

I 

I 



\ 



\ 

\ 



1 7.2 



1 7.4 
UTC (hrs) 



17.6 



17. 



Fig. 1. Bright Star Scan Results on ip— Ursa Major: The counting rates are shown 
as a function of time (UTC). The vertical lines around the peak of the profile are (i) 
expected transit time (solid line) (ii) centroid of the count rate profile (dot dashed 
line) defined by two dashed lines on either side. The baseline and the FWHM are 
also shown. 

Table 1 

Bright Star Scan Results of ijj— Ursa Major for a telescope 



Mirror 


Peak Count rate (Hz) 


FWHM (deg.) 


Offset (deg.) 


A 


9430 


2.29 


-0.23 


B 


1260 


2.46 


0.03 


C 


847 


2.50 


-0.14 


D 


781 


2.80 


0.09 


E 


493 


2.67 


-0.07 


F 


616 


2.18 


-0.19 


G 


2070 


2.38 


0.24 



This method of alignment is similar to that adapted for aligning the imaging 
telescope arrays [18]. 
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4 Estimation of Timing Resolution 

In order to estimate the angular resolution one has to estimate the error in 
timing measurement. To determine the arrival time of photons accurately, the 
PMT should have a fast response time and low dispersion in the photo-electron 
transit times. We use a low noise, high gain tube having rise-time of ~ 2 ns, 
pulse width of 3 ns and transit time jitter of 1.9 ns (FWHM). The PMT 
signals are brought to the recording stations through low loss coaxial cables 
(RG 213). The intrinsic timing jitter of the signals from the PMTs limit the 
accuracy of timing measurements. 

The timing resolution is estimated from the data. The relative delays of 6 
peripheral PMT's of each telescope and royal sum pulses with respect to the 
trigger pulse were recorded using the fast time to digital converters(TDC). The 
signals were delayed suitably using ECL based delay generators having an in- 
trinsic stability of ±0.1ns. Linearity of TDC was monitored by periodically 
calibrating them with standard cable delays. The LeCroy(Model # 2228A) 
and Philips Scientific (Mo del # 7186) TDC modules that are used in the ex- 
periment are set to a sensitivity of 0.25 and 0.2 ns per count, respectively. 

Data are collected with all telescopes pointing to zenith. The width of the 
distribution of differences in arrival times of signals (5Uj) of two TDC channels 
is an indication of the limiting accuracy of timing measurement, provided the 
signals originate from PMT's located nearby. The distribution of Stij tends to 
broaden due to one or more reasons mentioned below: 

(a) fluctuations in the arrival time and density of Cerenkov photons. 

(b) errors in mirror pointing result in sampling different angular regions, 
which in turn leads to pulse slewing. 

(c) spatial separation between the mirrors/telescopes. The mean arrival time 
as well as its fluctuations is a function of core distance [14]. 

(d) fluctuations in electron transit time in PMT's, the differences in the cable 
delays and the propagation delays in the electronic circuits. 

Delays due to unequal cable lengths, and propagation delays in the electronic 
circuits have been matched to an accuracy of about 0.25 ns. Mirrors are aligned 
to within ± 0°.2 as mentioned earlier. The effect of slewing responses of the 
PMT's is often determined by using filter wheels to attenuate laser calibration 
light flashes transported to each mirror through optical fibers [19]. No such 
measurements were made in this case since the pulse rise time is still reasonably 
small (< 2.4 ns). The effects due to finite separation between the detectors 
is corrected as discussed below and shown in figure 2. 

In the case of signals from individual PMT's of a telescope the distances 
between the detectors is of the order of a metre. Only those combinations 
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Table 2 

Timing Resolution of PACT Array 



Timing 
Information 


Sector 
# 


Timing 

Resolution(cjj ns) 


Royal Sum 


3 


2.3±0.1 


Royal Sum 


4 


2.1±0.2 


Royal Sum 


3+4 


2.2±0.1 


Individual PMT 


3 


l.liO.l 


Individual PMT 


4 


0.8±0.1 



corresponding to neighboring PMT's are considered in order to minimize the 
contributions from fluctuations in the arrival time of Cerenkov photons and 
distance dependent effects, to the uncertainty in timing measurement. The 
results are shown in the table 2. The limiting accuracy of timing measurement 
is about 1 ns for individual PMT signals. 

In the case of royal sum pulses, the distance between neighboring telescopes 
is 20 m or more, which is too large to be ignored. Moreover, the spread in the 
arrival time depends upon how well the signals of a telescope are added. To 
remove the distance dependent effects, the standard deviation are plotted 
as a function of separation between the telescopes, i and j and is shown in 
figure 2. In the past, the latter effect was minimized by tilting the telescopes 
towards each other by about a degree and the timing resolution was measured 
[7]. However we feel that the present method of reducing contribution from 
the detector separation is more accurate. From the figure it can be seen that 
the error aij clearly increases with the separation as expected. The FWHM 
for two telescopes side by side is estimated by extrapolation, thus removing 
the effect of finite separation between the telescopes on the distribution. 
The corresponding values are shown in table 2. These values are comparable 
to the intrinsic timing jitter for proton primaries estimated from simulation 
studies for PACT telescopes [15]. This means that the contribution to timing 
resolution from instrumental effects is negligibly small. 



5 Estimation of Arrival Direction of Shower 

The arrival direction of a shower is determined by reconstructing the shower 
front using the relative arrival time of Cerenkov photons at each telescopes. 
We fit the Cerenkov photon front to a plane using the measured arrival times, 
the normal to this plane gives the direction of shower axis. 
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Fig. 2. Figure shows the variation of < o~ij > /V2, derived from the TDC differ- 
ence distributions of royal sum pulses, as a function of detector separation for (a) 
single sector 3) and (b) for 2 sectors. Fitted straight lines extrapolated to zero 
separation yield the required timing resolution independent of telescope separation. 

5.1 Calculation of Time-offsets 



The relative arrival time of pulses as measured in the experiment is not the 
relative arrival time of Cerenkov photons at the PMT, needed for reconstruct- 
ing the shower front. A finite but constant delay between pulses from different 
PMT's arise due to unequal cable lengths, differences in electronic propagation 
delays and differences in photo-multiplier transit time etc. These are termed 
as TO or time-offsets. Thus the measured relative arrival times have to be cor- 
rected for this time-offsets to get the relative arrival time of Cerenkov photons 
at the PMT. In the absence of a calibration light source which could generate 
pulses in all PMT's simultaneously, these time-offsets have been determined 
from the data itself [20]. 

The time difference between Cerenkov photons arriving at PMT's (or tele- 
scopes) located side by side should be zero on the average for vertical showers. 
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Thus the average of time delays between two PMT's, from a large sample 
of data, is entirely due to difference between the two time-offsets. We chose 
vertical showers in order to eliminate the errors in the estimation of offsets 
due to finite telescope separations for inclined showers. 

If TOj and TOj are the time offsets for the PMTs i and j, we can write an 
equation of the form 

(TO, - TO,) = Qj (3) 

where CV, is the mean delay between a pair of PMT's i and j after correcting 
for the time difference due to differences in height (z-coordinates) of PMT's(or 
telescopes) . 

The total x 2 can De expressed as 

n 

X 2 = E WyiTOi-TOj-Cy) 2 (4) 

where n is the total number of PMT's, Wi/s are the statistical weight fac- 
tors, (Wij = l/(T C ij , where a c \j is the uncertainty in determining Cy). x 2 
minimization will give a set of n equations of the form 

n n 

E Wy(T0 i -T0 i )= E W t3 C t3 = C[ (5) 

3=^+3 3=l;i¥=3 

Thus we have a set of simultaneous equations involving W, TO and C. All 
these sets of equations can be written in the form of matrix equation 

W TO = C (6) 

Solving this matrix equation one gets an estimate of time-offsets for (n-1) 
detectors with respect to n th detector. It is to be noted that the differences in 
time delays due to the differences in the heights of the telescopes were taken 
off while calculating these time-offsets, as already mentioned. 

However the mean time difference < St > (ns) is a function of core distance 
R (m) and varies, for vertical showers, as 



< St >= 




where Sd is the separation between the telescopes, H the radius of curvature 
of the Cerenkov shower front and c the velocity of light. Typically, < St >= 
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0. 0117/? ns for the energy and altitude concerning us (H ~ 10 km). In our 
experiment, the average separation between the telescopes in a sector is 35 
m-50 m whereas the separation between the PMT's in a telescope is of the 
order of a metre. This can introduce a systematic error in the estimated values 
of T. For a proton shower of average energy (~ 1.8 TeV) the mean collection 
radius is ~ 100 m. Using this it is estimated that the maximum systematic 
error in T arising from the Cerenkov shower front curvature is ~ 1.2 ns. 

5.2 Reconstruction of Arrival Direction 

Using the plane front approximation the arrival direction of the shower is 
estimated as follows [20,10]. If Xi,yi,Z{ are the coordinates of the i th PMT, 

1, m,n the direction cosines of the shower axis and U the arrival time of the 
photons at this PMT then the equation relating them is given by 

Ixi + myi + nzi + c(U - t ) = (7) 

where to is the time at which the shower front passes through the origin of the 
coordinate system. Then the arrival direction of the shower can be estimated 
by a x 2 minimization where 

n 

X 2 = Y1 w S x i + m Vi + nz i + - ^o)) 2 (8) 
i=i 

where lOj's are the statistical weight factors for the i th timing measurement(= 
1/of where o~i is the uncertainty in the timing measurement for the i th detec- 
tor). The values of/, m, n and to are calculated using the equations dx 2 /dl = 0, 
dx 2 /dm = 0, dx 2 /dto = and I 2 + m 2 + n 2 = 1. We have ignored terms con- 
taining dn/dl and dn/dm as these are small. The contribution to the error 
in angle determination is negligibly small up to an angle of ~ 80° under this 
assumption. 

In the first iteration all PMT's with valid delays are used for the fit and the 
observed - expected delays are calculated for each. In the subsequent iteration, 
the PMT with the largest deviation is rejected if the absolute difference be- 
tween observed and expected delays is greater^ than 3 ns. This is to minimize 
the effect of isolated large fluctuations in the photon arrival times [21]. Fig- 
ure 3 shows the observed - expected delay distributions for various telescopes, 
based on royal sum pulses. This process of fitting is continued until there are 
no more large deviations or the number of available PMT's have become less 
than 4, in which case the shower front reconstruction process is aborted. 

1 ±3cr, The standard deviation of observed - expected delay ~ 1 ns 
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Figure 4 shows the zenith and azimuthal angle distributions of the recon- 
structed arrival directions using the procedure explained above. The recon- 
struction of the shower front was carried out using 25 telescope data collected 
with telescopes in the vertical position. The observed distributions are consis- 
tent with those expected. 






Fig. 3. Observed - Expected Delay Distributions for various telescopes in a sector, 
based on royal sum pulses. It is seen that the a of the delay distribution is about 1 
ns. 

Further, it has to be ascertained that the telescopes are pointing to the true 
zenith direction. A tilt in the observer's co-ordinate system will imply a sys- 
tematic error in pointing to the source direction. It is ensured that the assumed 
horizontal plane is really perpendicular to the direction of true zenith from 
the projected angle distributions of the fitted directions. The projected angles 
a and (5 are related to the zenith angle 9 and azimuth by 

tan a = tan 9 cos <fi 

and 

tan P = tan 9 sin <f> 



where a and (3 are the projected angles in the N-S and E-W plane respectively. 
The frequency distributions of a and (3 are shown in figure 5. The ratio of the 
number of events above and below zero are 1.0 ± 0.03 and 0.96 ± 0.03 for 
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Fig. 4. Distribution of zenith and azimuth angles estimated from royal sum timing 
informations from all telescopes. 

a and (3 respectively. The symmetry of these distributions thus verify the 
orthogonality of the telescope pointing direction to the horizontal. 



6 Angular Resolution of PACT 



We define the angular resolution of the array by the angle which corresponds to 
the 68% acceptance of mainly proton events. This definition is consistent with 
that for the HESS array except that their definition refers to 7-ray showers 
only [22]. The angular resolution of PACT has been estimated by using the 
split array method. The array is divided into two independent parts and the 
arrival direction is estimated for each shower from these two sub-arrays. The 
distribution of space angle (-0) between these two estimates is a measure of 
the angular resolution, i.e. the accuracy with which one can estimate the 
arrival direction. Since there are two independent estimates of the direction, 
the angular resolution is given by the peak[^] of the distribution of space angle 
between the two directions as ^S^. 



2 for small angle approximation of the point spread function coincides with the 
peak of the space angle (ip) distribution 
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Fig. 5. Distribution of Projected Angles, a and (3, in the (a) N-S and (b) E-W 
planes. The dot-dashed lines are Gaussian fits to the distribution 

Similarly, the uncertainties in the zenith angle 9 and azimuth <fi are estimated 
from the distributions of A9 and A0, the differences between the correspond- 
ing angles from two independent estimates for a shower. Figure 6 shows these 
distributions for A9 and sin6>.A0, respectively. 



6.1 Using Royal Sum pulses 



The array is divided into two independent parts of equal number of telescopes, 
say 12 telescopes from sectors 1,2 and sectors 3,4. The arrival direction is es- 
timated for each shower from these two independent arrays using the relative 
arrival times of royal sum pulses. The distribution of the space angle between 
the two estimated directions is shown in figure 7. Several combinations were 
tried to understand the dependence of angular resolution on the separation 
between the telescopes and/or the number of detectors. The results are sum- 
marized in table 3. The term 'odd-even' refers to 3 telescopes each from Sector 
3 and 4 grouped into one set of 6 telescopes and the remaining 6 grouped into 
another. 

It is seen from table 3 that the angular resolution improves as the separation 
between the detectors increases and also as the number of detectors increases, 
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Fig. 6. The differences of two independent estimates viz. AO and Acfi distributions 
using 6 telescopes for each estimate. 



Table 3 

Angular Resolution of PACT using Royal Sum pulses. 



No.of 
telescopes 


Combination of 
Detectors fitted 


Average 
detector 
separation (to) 


Peak of Space 
Angle Difference 
Distribution(deg. ) 


Angular 

Resolution 

(deg.) 


6 


Sector 3 vs 4 


31.8 


0.88 


0.62 


6 


Odd us Even 


48.8 


0.64 


0.45 


12 


Sector 3 and 4 vs 
Sector 1 and 2 


40.1 


0.59 


0.42 


12 


Sector 1 and 3 vs 
Sector 2 and 4 


52.0 


0.47 


0.33 



as expected. The dependence of angular resolution on D, the distance between 
telescopes and n, the number of telescopes used in the fit are shown in figure 8. 
The panel (a) of the figure shows the variation of angular resolution with D for 
(i) 4 telescopes and (ii) 12 telescopes (• and * respectively). The dependence 
on D is calculated to be ~ £>-o.94±o.o8 anci a ^ ^ D -i.o5±om reS p ec ti V ely, 

which is compliant with theoretically expected 1/D dependence. 
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Fig. 7. Distribution of space angle between two measured directions. The estimates 
are made with 12 telescopes from sector 1,2 and 12 from Sector 3,4. The dot-dashed 
line is a fit to the observed distribution. 

Once the 1/D dependence is established, the angular resolution is then cor- 
rected for the distance to extract the dependence on the number of degrees 
of freedom. This is shown in the panel (b) of figure 8. It is seen that ~ 
n -o.4i±o.n^ cons i s tent with what one expects, on a simple statistical basis, that 
cty will decrease as 

It is seen from table 3 that the best value of angular resolution, obtained using 
12 royal sum pulse information, is 0°.33. Since we have 25 telescopes, there 
will be a further improvement in the angular resolution by a factor of \/2 when 
the arrival direction is estimated using all the telescopes in the array. Thus 
the angular resolution of the array is estimated to be 0°.24, from the royal 
sum information of all telescopes. The corresponding threshold energy of the 
primary 7— ray is about 1.5 TeV. 



6.2 Using Individual PMT timing pulses 



The information on the relative arrival times of Cerenkov photons at the 
individual PMT's are available only within a sector. Only those events with 
information in all 6 telescopes in sectors 3 and/or 4 are used for estimating 
angular resolution. It is obtained by dividing the data from the sectors into 
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Fig. 8. Dependence of angular resolution a^, as a function of (a) separation between 
detectors and (b) the number of telescopes. 

2 subsets and obtaining the space angle between the two arrival direction 
estimates from the two subsets of data. The distributions of deviations of 
zenith and azimuth angle estimates are shown in figure 9 which corresponds 
to data from sector 3. 

Information for the 3 PMT's(labelled A, C and E) are grouped into one set 
while those for the remaining 3 PMT's of a telescope(i?, D and F) are grouped 
into another set. The results are summarized in table 4 . The rows 1 and 4 
correspond to the cases in which two subsets are obtained by demanding all 
A PMT's of 6 telescopes, all B PMT's etc. Rows 2 and 5 correspond to cases 
in which the 1st group consists of at least one valid TDC information in 
any of 3 A, C, E PMT' and the second group from any one of 3 B, D, F 
PMT's. Similarly the rows 3 and 6 correspond to two sets with at least 2 
valid TDCs in each telescope. Finally row 7 refers to the case in which the 
arrival directions are obtained separately from sector 3 and sector 4 events 
and collating event arrival times to pick common events. Column 2 shows the 
corresponding number of detectors used in the fit for all cases. Figure 10 shows 
the distribution of space angle between the two direction estimates made using 
individual PMT signals in a sector. 

A better estimate of arrival direction of the shower is possible with the use 
of individual PMT signals as seen in table 4. The angular resolution obtained 
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Fig. 9. The differences of two independent estimates viz. AO and smO.Acp distribu- 
tions using individual mirrors of sector 3. 

using individual PMT's of just one sector is 0°.24 (row 6 of table 4). When one 
considers the entire array the average separation doubles while the number of 
degrees of freedom quadruples. Consequently, there will be an improvement 
by a factor 4 in the accuracy of the arrival angle estimate. Based on these 
considerations, the angular resolution of the entire array is estimated to be, 
o$ ~ 0°.06., for > 2 active PMT's per telescope^ 

We have 6 timing information from each telescope. If we use the timing in- 
formation from a totality of 150 PMT's in the entire array, then the angular 
resolution works out to be: 

~ 0°.04 

The same figure is arrived at when we use row 7 of table 4. This estimate is 
based on the fact that > 4 PMTs have valid TDC information in a telescope. 
The primary energy threshold for such events is relatively higher, about 3 
TeV. The errors in the arrival direction estimate are larger for lower energy 
showers. While the angular resolution of the imaging telescopes are limited by 
the PMT sizes, only the future imaging telescope arrays (like VERITAS [23] or 
HESS [24]) could achieve a better angular resolution. PACT is able to achieve 

3 It may be noted that the average number of timing signals from individual mirrors 
in a sector is 15 when we demand > 2 PMT's and 27 when we demand > 4 PMT's 
in each telescope. 
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Table 4 

Angular Resolution of PACT using Individual PMT Information 



Sector 


No. of Detectors 


Combination of 


Peak of Space 


Angular 


# 


used in the fit 


Detectors 


Angle Distri- 


Resol- 








bution (deg) 


ution (deg.) 


3 


6 


all A, all B, etc 


0.46 


0.33 


3 


> 6 


at least 1 in 


0.45 


0.32 






each telescope 






3 


> 12 


~2 in a telescope 


0.43 


0.3 


4 


6 


all A, all B, etc 


0.48 


0.34 


4 


> 6 


at least 1 in 










each telescope 


0.39 


0.28 


4 


> 12 


~2 in a telescope 


0.34 


0.24 


3 and 4 


> 25 


~ 5 in a telescope 


0.2 


0.14 



this because it is a distributed array with multiple sampling of Cerenkov light 
pool. 



7 Discussions 

For a wavefront sampling experiments such as PACT, the technique by which 
off-axis events are rejected efficiently using their fine angular resolution, is 
vital. Hence significant effort has gone in to the improvement of the angular 
resolution of such arrays. In the case of the imaging telescopes the angular 
resolution is mainly decided by the size of the photo-tubes in the imaging 
camera at the focal point. However a new technique of stereoscopic viewing 
has been developed in 1996 by HEGRA group wherein the same shower is 
imaged by multiple telescopes. Using this technique the shower axis can be 
completely reconstructed in space using the information concerning the lo- 
cation of the images within each cameras and their angular orientation [18]. 
The reconstruction procedure allows determination of the shower direction, 
the core location and the height of the shower maximum. This technique has 
greatly improved the angular resolution of imaging telescope arrays at the 
cost of reduced collection area. On the other hand the improvement of angu- 
lar resolution of a wavefront sampling array results in a direct improvement 
in the flux sensitivity since the already increased collection area of such ar- 
rays is unaffected by the angle reconstruction. It may be noted here that a 
3-dimensional reconstruction of the shower is possible for the a non-imaging 
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Fig. 10. Distribution of space angle obtained from Sector 4 data, at least 2 PMT in 
a telescope have to have valid TDC 

array by making use of the curvature of the Cerenkov light front [25]. 

From the results of the present study it is seen that two independent mea- 
surements of the angular resolutions are ~ 0°.24 and 0°.06 resulting from the 
timing measurements of telescopes and mirrors of the entire array respectively 
for similar energy events. The essential difference between the two measure- 
ments is an increase in the number of degrees of freedom by a factor of ~ 2.5 
in the latter case. From this consideration the expected angular resolution for 
individual mirrors is 0.24/V2.5 = 0°.15 which is almost twice the value of 
0°.06 mentioned above. These are consistent since the timing uncertainty for 
royal sums is higher by the same factor. 

Figure 11 shows differential plots of the number of events per solid angle ^ 
as a function of the zenith angle. The panel (a) shows plots for events whose 
arrival directions have been estimated using timing signals from royal sums 
while the panel (b) shows the same for individual mirror timing information. 
The two curves shown in panel (a) correspond to events detected simultane- 
ously in all the 4 sectors (solid line) and in one sector (# 3) only (dashed line). 
This provides a method of measuring the effective opening angle of the array 
[27]. From these plots it can be seen that the array field of view is around 
2°. 3. This is significantly lower than the geometric opening angle of ~ 2°. 9. 
It may be noted that the effective field of view of the array is expected to be 
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Table 5 

Angular resolutions of various atmospheric Cerenkov experiments in the world. The 
asterisk indicates that the quoted angular resolution refers to this primary energy. 
In the rest of the cases the quoted angular resolutions are not necessarily applicable 
at the 7-ray threshold energies shown in column 4. In the case of imaging telescopes 
it may be noted that the actual angular resolution could be worse than the pixel 
size quoted here. 





Observatory 


Angular 


7-ray Energy 


Reference 




Name 


Resolution (deg.) 


Threshold (GeV) 




Imaging 


Whipple 


0.14 (€ 68%) 


500 


[28] 


Cerenkov 


CAT 


0.11 


250 


[29] 


Telescopes 


HEGRA(CTl) 


0.14 (€ 68%) 


1700 


[30] 




CANGAROO 


0.12 


1000 


[31] 




TACTIC 


0.31 


700 ± 200 


[32] 




SHALON 


0.4 


1000* 


[33] 




ALATOO I or II 








Stereoscopic 


HEGRA 


0.1 


500 


[34,35] 


System 


SHALON 
ALATOO 


0.1 


1000* 


[33] 


Non-Imaging 


THEMISTOCLE 


0.63 (€ 75%) 


3000 


[36] 


Cerenkov 


STACEE 


0.25 (€ 68%) 


190 ± 60 


[19] 


Telescope 


CELESTE 


0.2 


60 ± 20 


[37] 


Arrays 


GRAAL 


0.7 (€ 63%) 


250 ± 110 


[38] 




PACT 


2.4' (€ 68%) 


3000* 






AIROBICC 


0.29 ± 0.05 (€ 68%) 


30,000 


[39] 


EAS 


MILAGRO 


1.0 


1000 


[40] 




TIBET III 


0.87 (€ 50%) 


3000 


[41] 


Proposed 


VERITAS 


2.4'(€ 68%) 


1000* 


[23] 


Imaging 


MAGIC 


0.1, 0.2 (€ 68%) 


10 


[42] 


Arrays 


HESS 


0.1 (€ 68%) 


40 


[24] 




CANGAROO III 


0.1 


100 


[43] 
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Fig. 11. A differential plot of the number of events per solid angle as a function of 
zenith angle using estimates from (a) royal sum timing signals and (b) individual 
mirror timing signals. The solid line histogram in (a) corresponds to royal sum events 
detected simultaneously in all the 4 sectors while that in dotted line corresponds to 
royal sum events detected in sector 3 only. (See text for details.) 

lower than the geometric value primarily because the event trigger efficiency 
falls with increasing core distance and zenith angle. The other possible fac- 
tors which could potentially reduce the effective field of view like the residual 
alignment errors and non-ideal optical quality of the mirrors are relatively less 
significant. The effective FWHM of the array (~ 2°. 3) for which the mean tele- 
scope separation is ~ 70 m is lower than that for a sector (~ 2°. 6) for which 
the mean telescope separation is ~ 32 m supporting the above argument. 

The data used for estimating the angular resolution of PACT reported here 
consists of cosmic rays. It is well known that the protons exhibit a larger (by 
almost a factor of 4) intrinsic timing jitter relative to 7-ray primaries [15]. 
While the contribution of systematic effects to the timing resolution (which 
is estimated to be small) is independent of the primary species, the decreased 
intrinsic timing jitter is expected to improve the angular resolution for 7-ray 
primaries significantly. 

Table 5 summarizes the angular resolution of imaging as well as non-imaging 
atmospheric Cerenkov telescopes currently in operation and under construc- 
tion. It can be seen that the angular resolutions of most of the single imaging 
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telescopes are modest, in the range 0°.l — 0°.7. Wavefront sampling exper- 
iments using the modified solar concentrator arrays like the STACEE, CE- 
LESTE, GRAAL etc., will not be able to exploit their angle reconstruction 
technique to reject a significant cosmic ray background primarily because of 
their rather small field of view, in the range 0°.24 — 0°.7 (FWHM) which are 
comparable to their angular resolutions. PACT, on the other hand, has the best 
angular resolution among the non-imaging experiments which can be used to 
enhance the signal to noise ratio significantly. PACT has been able to achieve 
this because of multiple fast timing measurements at each telescope. This has 
two-fold advantages: firstly, it provides an increased number of degrees of free- 
dom which in turn improves the accuracy of angle reconstruction and secondly, 
it provides a means of computing the Cerenkov photon arrival time dispersion 
at each telescope. The timing jitter is a species sensitive parameter and hence 
will enable us reject a significant fraction of on-axis hadronic showers. The 
signal to noise ratio could also be improved by using the normalized x 2 values 
resulting from a spherical fit to the timing information of a shower. This is 
expected to be larger for hadron initiated showers [44]. 

The future imaging telescope arrays are expected to achieve unprecedented 
angular resolution by employing the stereoscopic technique. They are able to 
reconstruct the shower in 3-dimensions using the multiple images in 3 or more 
imaging telescopes. From table 5 it can be seen that both the HEGRA and 
the SHALON-ALATOO arrays are able to improve the angular resolution 
of their single imaging telescopes by a factor of ~ 4 by using 2 or more of 
them in stereoscopic mode. It may be mentioned here that the quoted angular 
resolutions for the future projects listed in the table are from simulation results 
for 7-ray primaries while the rest are derived from measurements from cosmic 
ray events. 

It is well known that the Cerenkov light front has a curvature. It was seen 
that when two well separated Cerenkov telescopes were tilted towards each 
other by about a degree the coincidence rate increased and also reduced the 
spread in the time separation between them [7,9]. This indicated, as claimed 
by the authors, the presence of curvature in the photon front. More recently, 
it was shown that the radius of curvature of the front is equal to the height 
of the shower maximum from the observation altitude [14]. Hence it is clear 
that a plane front approximation of the Cerenkov light front will introduce 
a systematic error in the arrival angle reconstruction. The large separation 
of telescopes, to some extent, offsets the worsening of resolution due to the 
curvature of the light front. It may also be pointed out that the effect of curva- 
ture on the angular resolution is more significant for near vertical showers. In 
the case of inclined showers, the arrival time differences at spatially separated 
telescopes due to shower axis inclination far exceed the differences arising out 
of the shower front curvature. However it has been argued that the angular res- 
olution of given array will improve if one corrects for the wavefront curvature 
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[26] . The details of the systematic effects due to the plane front approximation 
and the improvement in the accuracy of estimated arrival directions when the 
curvature of the shower front is taken into account are currently under in- 
vestigation. A paper based on the results of these simulation studies is under 
preparation. 



8 Summary 

A detailed analysis of the angular resolution of PACT using data collected with 
telescopes pointing to zenith is presented. The improvement in the angular 
resolution with larger separation between detectors (D) and with increase in 
the number of degrees of freedom (n) has been verified. The angular resolution 
cty, is found to improve according to the relation: 

1 

D^Jn 

There are two types of angular resolutions that could be defined for PACT. 
Firstly, the fast timing information from the 25 telescopes spread in an area 
of about 80 m x 100 m yield an angular resolution of ~ 0°.2. Secondly the 
timing information from the mirrors constituting the 25 telescopes yield the 
best angular resolution of ~ 2.4' at a 7-ray energy of around 3 TeV. This is the 
best angular resolution achieved for any ground based atmospheric telescope 
system which will probably be superseded by the proposed VERITAS or HESS 
imaging telescope array. 

The angular resolution for 7-ray primaries could be significantly better than 
what is presented here, which is based on cosmic ray primaries, because of 
two main reasons. Firstly, photon arrival time jitter for 7-ray primaries is far 
less than that for charged cosmic rays and secondly the radius of curvature for 
7-ray primaries is more than that for cosmic rays which consequently reduces 
the effects of the shower front curvature for the former. 
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